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Magnetic helix wire is one of the most simple magnetic systems which manifest properties of both 
curvature and torsion. There exist two equilibrium states in the helix wire with easy-tangential 
anisotropy: a quasi-tangential magnetization distribution in case of relatively small curvatures and 
torsions, and an onion state in opposite case. In the last case the magnetization is close to tangential 
one, deviations are caused by the torsion and curvature. Possible equilibrium magnetization states 
in the helix magnet with different anisotropy directions are studied theoretically. The torsion also 
essentially influences the spin-wave dynamics, acting as an effective magnetic field. Originated 
from the curvature induced effective Dzyaloshinskii interaction, this magnetic field leads to the 
coupling between the helix chirality and the magnetochirality, it breaks mirror symmetry in spin- 
wave spectrum. All analytical predictions on magnetization statics an dynamics are well confirmed 
by the direct spin lattice simulations. 

PACS numbers: 75.30.Et, 75.75.-c, 75.78.-n 


I. INTRODUCTION 

During the past few years there is a growing inter¬ 
est to curvature effects in physics of nanomagnetism. 
A crucial aspect of the interest is caused by recent 
achievements in nanotechnologies of flexible, stretchable 
and printable magnetoelectronics (see Ref. 1 and refer¬ 
ences therein). Effects of the curvature on the magne¬ 
tization structure in nanomagnetic particles of nontriv¬ 
ial geometry were studied for cylinders,torus,half- 
spheres^, spherical shells,^ hemispherical caps,^’^ cylin¬ 
drical capped nano membranes,^, cone shells^^’^^, and 
paraboloidal shells^^. Chiral and curvature effects with 
account of the nonlocal dipolar interaction were discussed 
for cylinder nanotubes. 

Very recently we have developed fully three dimen¬ 
sional (3D) approach for studying statics and dynamics 
of thin magnetic shells and wires of arbitrary shape. 

This approach gives a possibility to derive the energy for 
arbitrary curves and surfaces and arbitrary magnetiza¬ 
tion vector fields on the assumption that the anisotropy 
contribution greatly exceeds the dipolar and other weak 
interactions, i.e. for hard magnets. We have shown 
that due to the curvature two additional effective mag¬ 
netic interactions originate from the exchange term: (i) 
curvature induced effective anisotropy, which is bilinear 
with respect to curvature and torsion and (ii) curvature 
induced effective Dzyaloshinskii interaction, which is lin¬ 
ear with respect to curvature and torsion. This novel 
approach open doors for studying several perspective di¬ 
rections in nanomagnets, including topologically induced 
patterns^’^^ and magnetochiral effects^^’^^. 

The simplest system which displays both the prop¬ 
erties of the curvature and torsion is a helix wire, 
which is characterized by constant curvature and tor¬ 
sion. The interest to such a geometry is motivated by 
recent experiments on rolled-up ferromagnetic microhe¬ 


lix coils.Depending on the anisotropy direction dif¬ 
ferent artificial complex helimagnetic-like configurations 
were experimentally realized: hollow-bar-, corkscrew-, 
and radial-magnetized 3D micro-helix coils.Rolled 
magnetic structures are now widely discussed in the 
context of possible application in flexible and stretch¬ 
able magnetoelectronic devices,in particular, rolled- 
up GMR sensors,for magnetofluidic applications, spin- 
wave filters,and microrobots^^. Helix coil magnetic 
structures have the potential to be used a variety of 
bioapplication areas, such as in medical procedures, cell 
biology, or lab-on-a-chip.^^ 

In the current study we apply our theory^ ^ aimed to 
describe magnetization statics and linear dynamics in the 
helix wire. We analyze equilibrium states for different 
types of magnetocrystalline anisotropy. The equilibrium 
state is determined by the relationship between the cur¬ 
vature, torsion and the anisotropy strength: we describe 
possible magnetization distributions analytically. For 
three types of anisotropy (easy-tangential, easy-normal 
and easy-binormal) we compute phase diagrams of pos¬ 
sible equilibrium states. In each of these cases the equi¬ 
librium state is either onion one (high curvatures) or 
anisotropy-aligned state (for small ones), these results 
are summarized in Fig. 4. For example, in the most 
interesting case of easy-tangential anisotropy a quasi- 
tangential magnetization distribution appears for strong 
enough anisotropy, see Fig. 3(a,b). We show that pure 
tangential magnetization distribution is impossible. The 
deviation from the tangential state is determined by the 
the curvature and torsion; besides there exists the cou¬ 
pling between the helix chirality and the magneto chiral¬ 
ity of magnetization distribution. 

We study the problem of spin wave dynamics in the 
helix wire. Our analysis shows that the curvature and 
torsion act on magnons in two ways: besides the stan¬ 
dard potential scattering of magnons, there appears an 
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effective torsion induced magnetic field. The vector po¬ 
tential of effective field is mainly determined by the prod¬ 
uct of the torsion and the magneto chirality. The origin 
of this field is the curvature induced effective Dzyaloshin- 
skii interaction.^^ Finally, the torsion breaks the symme¬ 
try of spin wave spectrum with respect to the direction 
of spin wave propagation, see Fig. 5 . This effect is com¬ 
pletely analogous to the effect of asymmetry of magnon 
dispersion due to the natural Dzyaloshinskii interaction 
in magnetic films. 

The paper is organized as follows. In Sec. II we intro¬ 
duce the model of the curved wire and discuss different 
anisotropy-aligned states. The model of the helix wire 
appears in Sec. III. Equilibrium magnetization distribu¬ 
tions are describe analytically for the easy-tangential he¬ 
lix wire: the quasi-tangential state (see Sec. Ill A) and 
the onion one (see Sec. IIIB). The phase diagram of en¬ 
ergetically preferable states appears in Sec. IIIC. The 
problem of spin-wave dynamics is discussed in Sec. IV. 
In Sec. V we study statics and linear dynamics for he¬ 
lix wires with other anisotropy orientations: the easy- 
normal anisotropy (Sec. V A) and easy-binormal one 
(Sec. V B). We verify our theory by numerical simulations 
of the helix-shaped chain of discrete magnetic moments 
in Sec. VI. In Section VII we present final remarks and 
discuss possible perspectives and generalizations, in par¬ 
ticular, how to take into account magnetostatics effects. 
Some details about the computation of the onion state 
are presented in Appendix A. 


II. THE MODEL OF A CURVED WIRE 


where u and v are coordinates within the wire cross sec¬ 
tion {\u\, |i;| < h). 

Let us describe the magnetic properties of the wire. 
The magnetic energy of the wire can collect different 
contributions such as energies of exchange interaction, 
anisotropy, and dipolar one. We start our analysis with 
the case of a hard magnet where the anisotropy contri¬ 
bution greatly exceeds the dipolar and other weak inter¬ 
actions. For such hard magnets a quality fact or 

Q=^ 

^ ~ 27 rM 2 

is supposed to be large; here A > 0 is the constant of 
magnetocrystalline anisotropy and Mg is the saturation 
magnetization. 

We assume the magnetization spatial one- 
dimensionality, which can be formalized as m = m(s, t). 
This assumption is appropriate for the cases when 
the thickness h does not exceed the characteristic 
magnetic length w = ^^AfK with A being the exchange 
constant. The wire thickness is also supposed to be 
small in comparison with radii of curvature and torsion. 
Therefore our model provides an adequate picture under 
the following assumptions: 

h<re<C-, -, Q > 1 . ( 5 ) 

K, r 

That is why in the current study one can restrict our¬ 
self to the consideration of Heisenberg magnets with the 
energy 


We consider a model of a curved cylindrical wire. Let 
7(s) be a ID curve embedded in 3 D space with s 
being the arc length coordinate. It is convenient to use 
Frenet-Serret reference frame with basic vectors Cq,: 




= 7 



63 — Gt X Cn (1) 


with Ct being the tangent, being the normal, and 
Cb being binormal to 7. Here and below the prime 
denotes the derivative with respect to the arc length s 
and Greek indices a, [3 numerate curvilinear coordinates 
(TNB-coordinate) and curvilinear components of vector 
fields. The relation between and is determined by 
Frenet-Serret formulas: 
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Here k and r are the curvature and torsion of the wire, 
respectively. 

The wire of a finite thickness h can be defined as the 
following space domain 


( 3 ) 


E — AS f ds(^ex + ^an )5 

^ 2 (6) 
^ —2 tp ’ ^an) 

(^x — ‘ ^ TfX^ ^an — n 

with Can being the unit vector along the anisotropy axis, 
and S being the cross-section area. 

Typically, orientation of the anisotropy axis Can is de¬ 
termined by the wire geometry, e.g. it can be tangential 
to the wire,^^ which means in general complicated spatial 
dependence due to the curvilinear geometry. Therefore 
it is convenient to represent the energy of the magnet 
in the curvilinear reference frame (1), where (^an has a 
simplest form. For an arbitrary thin wire the exchange 
energy density can be presented as follows^ ^ 

^ex = ^ex + Cx + ^ex5 Cx=l'^^l ^ 

4x = Kapnianip, = Fap {niam’i^ - m'^mp) . 

Here the first term describes the common isotropic 
part of exchange expression which has the same form as 
for the straight wire. The second term describes an 
effective anisotropy interaction, where the components of 
the tensor Kai3 = FauE^u are bilinear with respect to the 
curvature k and the torsion r. This term is similar to the 
“geometrical potential” Note that a curvature caused 


r{s, u, v) = 7(s) + ue^ -h 
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Anisotropy 

type 

Anisotropy 

axis 

Gan 

Magnetization states 
in a helix wire 

Equilibrium states Orientation 
according 
to Ref. 17 

Easy-tangential 

Gt 

quasi-tangential corkscrew 
and onion 

Easy-normal 

Gn 

normal and onion radial 

Easy-binormal 

Gb 

quasi-binormal hollow-bar 

and onion 


TABLE L Types of equilibrium magnetization states for var¬ 
ious uniaxial anisotropies in a helix-shaped magnetic wire. 


“geometric” effective magnetic field was considered re¬ 
cently for curved magnonic waveguides.The last term 
in the exchange energy functional is the curvature 
induced effective Dzyaloshinskii interaction, which is lin¬ 
ear with respect to curvature and torsion. We will see 
below that this effective interaction causes an effective 
magnetic field; namely this interaction is responsible for 
the magnetochiral effects. 

We consider three types of curvilinear uniaxial 
anisotropy which correspond to possible curvilinear di¬ 
rections (1), see Table I: (i) an easy-tangential anisotropy 
corresponds to the anisotropy axis Can directed along Ct, 
where the anisotropy interaction tries to orient the mag¬ 
netization along the curve. Note that in soft magnets 
such kind of anisotropy appears effectively as a shape 
anisotropy caused by the dipolar interaction.^^ (ii) An 
easy-normal anisotropy is determined by the normal vec¬ 
tor Cn- (iii) An easy-binormal anisotropy direction cor¬ 
responds to the binormal basic vector e^. 

All three types of anisotropic magnets can be real¬ 
ized experimentally: In straight nanostrips/nanowires 
the anisotropy can have well-defined uniaxial directions, 
e.g., in-plane along the strip, in-plane perpendicularly 
to the strip, or out-of-plane, which corresponds to the 
uniformly magnetized samples in the corresponding di¬ 
rection. Using the coiling process,it is possible to 
obtain 3D microhelix coil strips with different magne¬ 
tization orientation: corkscrew-, radial-, and hollow-bar- 
magnetized, see Table I to get a link between anisotropy 
type and the magnetization orientation. 

For the further analysis it is convenient to introduce 
the angular parametrization of the magnetization unit 
vector m using the local Frenet-Serret reference frame: 

m = sin 6 > cos 0 Ct + sin 0 sin 0 Cn + cos 0 Cb, ( 8 ) 


where angular variables 0 and (j) depend on both spatial 
and temporal coordinates. Then the energy density (7) 
reads: 


^ex= — Tsin0]^+ [sin6>(0' k) — r cos6>cos0]^ (9a) 

sin^ 0 cos^ <p sin^ 0 sin^ 0 
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FIG. 1. (Color online) Schematics of the helix wire of the 
radius R and the pitch P. (a) Arrangement of the curvilinear 
Frenet-Serret reference frame (gt, Gn, Gb) from the front view, 
(b) Arrangement of the magnetization angles 0 and cj) with 
respect to the magnetization unit vector m. 


Here and denotes anisotropy energy 

densities of easy-tangential, easy-normal, and easy- 
binormal types, respectively. 

The magnetization dynamics follows the Landau- 
Lifshitz equation. In terms of angular variables 0 and 
(j) these equations read 


—^ sin Odtcj) 
7o 


5E 

Jo' 


-^ sin OdtO 

7o 


SE 

6(f) 


( 10 ) 


with Ms being the saturation magnetization and 70 being 
the gyromagnetic ratio. 


III. EQUILIBRIUM MAGNETIZATION STATES 
OF A HELIX WIRE WITH EASY-TANGENTIAL 
ANISOTROPY 

We study the curvilinear effects using the helix geom¬ 
etry, which is the simplest geometry which manifests the 
properties of both curvature and torsion. A typical pa¬ 
rameterization of the helix wire reads 

7 (x) = xRcosxEyRsinxE zpx, (Ha) 

where R is the helix radius, p = Pj (27r) with P being the 
pitch of the helix, and x is azimuthal angle of a cylindrical 
frame of reference with £-axis aligned along the helix 
axis, see Fig. 1. The helix has the constant curvature 
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K = R/{R^ and the torsion r = p/{R^ +P^)- For 

the further analysis it is instructive to rewrite ( 11 a) as 
a function of the arc length s and in terms of curvature 
and torsion 


7 (s) = xns^ cos ( — ) + sin ( — 

1 

One has to notice a one-to-one correspondence between 
(i?,p)-parametrization ( 11 a) and {K.,r)-one ( 11 b). 

In order to derive the explicit form of Landau-Lifshitz 
equations, we substitute the energy functional (9) into 
the Landau-Lifshitz equations (10): 


+ ZSqTS, 


( 11 b) 


M, 


sin Odtcj) = T cos (j) (/^ cos 20 — 29s(/) sin^ O) 

iS^an 


270^1 
+ dssO — sin 6 >cos 6 > 
M< 


(/^ + dgcj)) — r^cos^ (j) 


2 dO ’ 

sin OdfO = sin 0 cos 0 [2ds0 {k, + 9s0) — kt sin 0 ] 

iS^an 


270^1 

+ sin^ 6 > [dss(l^ + 2rds0 cos (p 


r^smpcoscj)] — -- 


( 12 ) 


where S’an is the density of the anisotropy energy, see 
(9b). 

We are mostly interested in the case of easy-tangential 
anisotropy, which is typical for the wires. In this case the 
anisotropy energy density has the form see (9b). 

First we discuss the limit case r = 0 (ring wire instead 
of the helix). For any plane curve the energy functional 
(9) with easy-tangential or easy-normal anisotropy is 
minimized by the plane magnetization distribution, Oq = 
7 r/ 2 . The energy minimization in respect to p results in 
the pendulum equation 

— sin^cos 0 = 0, >c = nw (13) 



FIG. 2. (Color online) Equilibrium magnetization distri¬ 
bution in the quasi-tangential state of the helix wire with 
C = +1. Lines correspond to the analytics, see Eq. (17). Sym¬ 
bols correspond to simulations: (a) anisotropic Heisenberg 
magnets [see (50)] (Q — 2^ w — 1)^ (am) wires with account 
of dipolar interaction [see (49), (53)] (Q = 2, = 2-£/\/5), 

and (ms) isotropic wires with account of dipolar interaction 
[see (49),(53)] (Q = 2, = 21). The boundary curve 0l{a) 

corresponds to (23). 


A. Quasi-tangential state 

Now we consider the helix wire with a finite torsion, 
r 7 ^ 0. Similar to the case of a ring wire, discussed above, 
we look for the homogeneous (in the curvilinear reference 
frame) solution. Such kind of solutions is possible due 
to the constant curvature k and the torsion r. We can 
easily solve the static equations, see Eq. (12), using the 
substitution 0{s) = 0^ and 0 (s) = 


with X being the reduced curvature. 

The equilibrium magnetization state of a ring is a ho¬ 
mogeneous (in the curvilinear reference frame) vortex 
state in case of relatively small reduced curvature 
K < xq ~ 0.657 and inhomogeneous onion solution 0°^ 
for X > 

0 vor_o, 7 r, 0°^ = ^ — am(x,/c), X = — K(/c). (14) 

2 TT 

Here am(x, k) is the Jacobi amplitude^^ and the modulus 
k is determined by condition 


tan 26>* = - -- ^ ^ , 7 = 0, tt, (16) 

where C = = ± 1 , the quantity a = wr is the 

reduced torsion. Explicitly for magnetization angles we 
get 

2 Ccrx . 

0 = — — arctan ———, <p = 0 , tt, 

2 Vo 

Vo = l + (j^ - + Vu 

Vi = y/(l — + cr^)^ + 4x^(7^. 


2KkK{k) = TT (15) 

with K{k) being the complete elliptic integral of the first 
kind.^^ 


The dependence 0^{k^(j) is presented in Fig. 2. 

In the limit case of very small curvature and torsion 
(x, cr ^ 1), the magnetization distribution becomes al¬ 
most tangential, see Fig. 3(a) with the asymptotic be- 
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O^^JLX. 


(17) reads 


S'* = - 


1 - x2 - 0-2 + Vi 

2w^ 


(19) 


It should be noted that the magnetization state in the 
helix nanowire was recently studiedin particular, the 
magnon spectrum was shown to be affected by the curva¬ 
ture, which acts mainly as effective anisotropy. However 
the equilibrium state was forcedly supposed to be the 
tangential one in Ref. 33. 


B. Onion state 


(b) Quasi-tangential state: x = cr = xo ~ 0.657 



(c) Quasi-tangential state: x = a = 1 



(d) Onion state: x = 1.5, cr = 1 

FIG. 3. (Color online) Magnetization distributions in the he¬ 
lix wire with C = -hi and easy-tangential anisotropy according 
to simulations data, see Sec. VI A. 


havior 

0^ ^ ^ for X, cr ^ 1. (18) 

That is why we refer to the state (17) as to the quasi- 
tangential state. Such a state is an analogue of the vortex 
state for the case of the torsion presence. 

Even in the strong anisotropic case the magnetization 
deviates from the tangential distribution: the inclination 
angle depends on the sign of Ccr. One can interpret the 
sign of cr as the helix chirality (different for right handed 
helix when a > 0 and left-handed one when a < 0); 
the quantity G can be interpreted as the magnetochiral¬ 
ity, hence on can say about coupling between the two 
chiralities. 

The energy density (9) of the quasi-tangential state 


Let us discuss the case of a large curvature and tor¬ 
sion. In analogy with the ring wire, we are looking for a 
solution periodic with respect to y, which is an analogue 
of the onion solution (14). Hence we look for solutions of 
the following form 

= f + ^(x)> = -X + (20a) 

with 'd(x) and (f(x) being 27r-periodic functions. Using 
an analogy with the ring case (a = 0) with exact onion 
solution (14) we name (20a) an onion solution. 

Numerically we found onion solutions for x > xq ~ 
0.657 in a wide range of cr, see Figs. 3(c), 4(a). The 
symmetry of the static form of Eqs. (12) dictates the 
symmetry of 27r-periodic functions 'd and (^, which has 
the following Fourier expansion 

N N 

= E cos(2n - l)x, (fix) = E 

n=l n=l 

(20b) 

where N ^ oc. By substituting series (20b) into the 
static version of Eqs. (12), one get the set of nonlinear 
equations for amplitudesand see (A5). Finally, 
the energy of the onion state <f°^(cr, x^), averaged over 
the helix period, can be calculated numerically using am¬ 
plitudes 'dn and ipn^ see Appendix A for details. 

C. Phase diagram 

Now we summarize results on the equilibrium mag¬ 
netization distribution. By comparing energies of differ¬ 
ent states, we compute the energetically preferable states 
for different curvature and torsion values. The result¬ 
ing phase diagram is presented in Fig. 4(a). There are 
two phases: (i) The quasi-tangential state is realized for 
relatively small curvatures, when x < X 5 (cr); in such 
a state the magnetization direction is close to the di¬ 
rection of easy-tangential anisotropy Ct, see Fig. 3(a,b) 
with the limit vortex orientation in case of the ring wire 
(r = 0). (ii) The onion state corresponds to the case, 
when X > X 5 (cr); the magnetization distribution is inho¬ 
mogeneous in accordance to (20), see Fig. 3(c). 
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EIG. 4. (Color online) Phase diagram of equilibrium magnetization states for the helix wire with different types of anisotropy. 
Symbols correspond to simulation data: green diamonds to homogeneous (in curvilinear reference frame) states and open circles 
to the onion ones, (a) Easy-tangential case, the curve Xb(cr) (solid green line), calculated by (21) with N = 3, describes the 
boundary between the quasi-tangential and the onion states; the dashed-dot line corresponds to ^^(cr) with N = 1. The curve 
Xc(cr) (dashed red line) describes the boundary of linear instability of the quasi-tangential state, the dotted line is its fitting 
by (33). In the region between lines ^^(cr) and Xc((7) the quasi-tangential state is metastable, (b) and (c) correspond to 
easy-normal and easy-binormal anisotropy, respectively; all notations have the same sense as in (a). Note that Eig. (b) has 
different scale in order to show the normal state region in details. 


The boundary between two phases = ^^(cr) can be 
derived using the condition 

( 21 ) 

where is energy density of the onion state averaged 
over the helix period 27rso, see (A6). The onion solution 
(20) is energetically preferable when its energy is lower 
than the energy of the quasi-tangential state (19). We 
computed the boundary curve numerically for = 1 and 
N = 3, see dot-dashed and solid lines, respectively in the 
Fig. 4(a). The obtained curves are very close, so the 
approximation = 1 is reasonable. This is because the 
onion state of the helix wire is very close to an uniform 
magnetization, see Fig. 3(c). 

For the approximate description of the boundary de¬ 
pendence we use the trial function 

+20-2, (22) 

which fits the numerically calculated curve ^^(cr) with an 
accuracy of about 5 x 10“^. 

Using the boundary dependence ^^(cr), one can easily 
compute domain of applicability of the quasi-tangential 
solution (17): 


IV. SPIN WAVE SPECTRUM IN A HELIX 
WIRE WITH EASY-TANGENTIAL 
ANISOTROPY 


We limit our consideration of spin waves by the case 
of the quasi-tangential magnetization state. First we lin¬ 
earize the Landau-Lifshitz equations (12) on the back¬ 
ground of the quasi-tangential equilibrium state (17), 

0{s,t) =e^ + (p{s, t) =(p* + ^ . (24) 

smt/^ 

Then for and (p we get the set of linear equations: 


dt>(p = + Vii9 — 2Ad^(p, 

= -d^^ip + V2ip + 2Ad^'d, 


where dt' is the derivative with respect to dimensionless 
time t' = Clot with CIq = 2K^/Ms and is the derivative 
with respect to dimensionless coordinate ^ = s/w. Here 
Vi is determined according to (17), the quantities V 2 and 
A have the following form: 


U 2 


1 + + Vi 

2 


A = -xcos6>^ - crQsmO^ 



(26) 


0^ e 


[$1, 1 ) , when Ccr > 0, 
(|, , when Ccr < 0. 


(23a) 


Here 0^ = 01(a) determines the boundary curve. 


While Vi and V 2 appear as scalar potentials, A acts 
as a vector potential A = Ae^ of effective magnetic 
field. This becomes obvious if we combine the set of 
linearized equations for and (p in a single equation for 
the complex-valued function ip = 1 } icp^ 


(23b) -idt'ip = H'tp + Wxp*, H = {-id^ - Af + U. (27a) 


6»b(cr) = e\>Cb{a),a). 
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This differential equation has a form of generalized 
Scrodinger equation, originally proposed for the descrip¬ 
tion of spin waves on the magnetic vortex background. 
The “potentials” in Eq. (27a) read 

2 ’ 2 2 ■ 

(27b) 

An effective magnetic field A is originated from the cur¬ 
vature induced effective Dzyaloshinskii interaction, see 
Eq. (7): the energy density harmonized using (24), 
reads^^ 

argt/). (28) 

Now we apply the traveling wave Ansatz for the spin- 
wave complex magnon amplitude 

'0(^, t') = ue*^ + 4> = + 77, (29) 

with q = kw being the dimensionless wave number, 
ft = uj/fto the dimensionless frequency, 77 is arbitrary 
phase, and u, v G M being constants. The corresponding 
wave vector is oriented along the wire, q = qe^; its ori¬ 
entation with respect to the equilibrium magnetization 
is determined by Eq. (17). By substituting the Ansatz 
(29) into the generalized Scrodinger equation (27), one 
can derive the spectrum of the spin waves: 

n{q) = 2Aq + ^(q^ + V,) {q^ + V 2 ). (30) 

Similar to the straight wire case with fls(g) = 1 + 
the spectrum of spin waves in the helix wire has a gap, 
caused, first of all, by the anisotropy (in dimensional 
units the gap has an order of IIq K). However its 
value essentially depends on the curvature and the tor¬ 
sion. Moreover, the spectrum gap occurs at finite g = go, 
see Eig. 5. This means the asymmetry in the spectrum 
with respect to the change g ^ —g: spin waves have 
different velocities depending on the direction (along the 
helix axis or in opposite direction). This asymmetry in 
the dispersion law (30) occurs in the first term 2Ag, which 
is originated from the effective Dzyaloshinskii interaction 

^ex • 

In this context it is instructive to mention that the 
spin wave spectrum in the presence of Dzyaloshinskii- 
Moriya interaction is known to be asymmetric with re¬ 
spect to wave vector inversion and has the minimum at 
finite wave vectors.The curvature induced asymme¬ 
try in the spin waves propagation in nanotubes and its 
analogy with the Dzyaloshinskii-Moriya interaction was 
discussed recently in Ref. 36. The spin-wave spectrum 
for the helix wire was calculated recently in Ref. 33, how¬ 
ever the deviations from the pure tangential state were 
no taken into account and the effective Dzyaloshinskii 
was not considered. 

In order to make analytical estimations, we consider 
now the dispersion law in case of very small curvatures 


and torsions: 

D(g) = (2gap + (g - Ccr)^ + o (x^, cr^, xcr) , 

^2 (31) 

One can see that the spin wave spectrum becomes asym¬ 
metrical one with increasing the curvature and the tor¬ 
sion: the minimum of the frequency corresponds to 
go = crC (in dimensional units the corresponding wave 
number = rC), its sign is determined by the product 
of the helix chirality and the magnetochirality. 

The further increase of the curvature and torsion de¬ 
crease the gap 12gap; there is a critical curve Xc = Xc(cr), 
where the gap vanishes, 0(gc) = 0 and dqQ{qc) = 0. 
One can easily find that qc = G\/A^ — U and the critical 
curve Xc = Xc(cr) can be found as a solution of algebraic 
equation 

AA^U = W‘^. (32) 

The critical curve Xc(cr), calculated numerically is plot¬ 
ted in Eig. 4(a) (dashed red curve). Eor the approximate 
description of the critical dependence we use the trial 
function 

= Vl + 2a2, (33) 

which fits the numerical results with an accuracy of about 
5 X 10“^, see the dotted curve in Eig. 4(a). In the region 
between the boundary curve and the instability 

curve Xc(cr) [see Eig. 4] the quasi-tangential state be¬ 
comes metastable. 


V. HELIX WITH OTHER ANISOTROPY 
ORIENTATIONS 

Let us discuss other types of anisotropies: easy-normal 
and easy-binormal, see Eq. (9b) and Table I. 


A. Easy—normal anisotropy 

Let us start the analysis of the easy-normal anisotropy 
with the limit case of the ring (r = 0). In this case, 
similarly to the easy-tangential anisotropy, the magneti¬ 
zation lies within the ring plane: 0 = 7r/2. The energy 
minimization with respect to (j) results in the pendulum 
equation: 

x^9^^0-h sin0cos0 = 0. (34) 

In analogy with the easy-tangential anisotropy the equi¬ 
librium state is the exactly normal state (j)'^ = ±7r/2 in 
case of relatively small reduced curvature x < xq and 
inhomogeneous onion solution 0^^(x) = 7r/2 — 0°^(x) for 
X > Xq, where function 0°^(x) is defined by (14). 










1 2 
Reduced wave number q 

(a) Easy-tangential anisotropy 


X = 0, (7 = 0 
— x = 0.25, o-e = 0.25 
■ - X = 0.45, o-e = 0.45 



0.5 1 

Reduced wave number q 

(b) Easy-normal anisotropy 


1 2 
Reduced wave number q 

(c) Easy-binormal anisotropy 


EIG. 5. (Color online) Top row demonstrates dispersion laws for spin waves in the helix wire for different anisotropies. The 
equilibrium states are homogeneous in the curvilinear reference frame. Symbols correspond to simulation data, see Sec. VIB, 
and lines to the analytics, see Eq. (30) and (39). Eew examples of dispersion relation are shown at the bottom row in terms of 
density plots to demonstrate that (30) is a single frequency branch in the system. 


In case of finite torsion there also exists exactly normal 
state 


0 ^ = 


TT 

2 ’ 





1 - 




(35) 


where C = ±1, see Fig. 6(a). Such a state is energet¬ 
ically preferable for relatively small values of >c and a. 
The magnetization in the normal state is directed ex¬ 
actly radially, which is well pronounced in experiments 
with 3D microhelix coil strips. 

In case of large curvature, there is the periodic (in 
curvilinear reference frame) onion solution, which has the 
form (20), see Fig. 6(b). Using the same numerical pro¬ 
cedure as in Sec. IIIB, we evaluate the onion solution 
and compute the phase diagram, see Fig. 4(b). 

For the approximate description of the boundary 
between two phases we use the fitting function 





0.67, 


(36) 




FIG. 6. (Golor online) Magnetization distribution in the helix 
wire with C = +1 and easy-normal anisotropy according to 
simulations data, see Sec. VIA. 


which fits the numerically calculated curve x^^(cr) with Let us discuss now the linear excitations on the back- 
an accuracy of about 3 x 10“^. ground of the normal solution. Using the same approach 
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as in Sec. IV, we linearize Landau-Lifshitz equations (12) 
on the background of the normal solution (35), 0 = 

(j) = (jf^ After linearization one gets a general¬ 

ized Scrodinger-like equation for the complex variable 
t/; = 'd + i(^, 

- idt^^ = (-% + V^) -h (37a) 

Here the “potentials” read 

-h 1 9 

C/" = l- W^ = -{e<j-ixf. (37b) 

Let US compare this equations with the generalized 
Scrodinger-like equation (27). First of all, there is no 
effective vector potential, since there is no asymme¬ 
try by effective Dzyaloshinskii interaction like in easy- 
tangential case. The second difference is that the poten¬ 
tial W in (37b) is a complex-valued one, hence the scat¬ 
tering problem is similar to the two-channel scattering 
process. Similar to (29) we apply the following traveling 
wave Ansatz for the spin-wave complex magnon ampli¬ 
tude 

iPi^, t') = tAie** + ^ = qi-nt' + rj, € C. 

(38) 

The difference is that constants are complex ones. 
Now by substituting the Ansatz (38) into the generalized 
Scrodinger equation (37), one can derive the spectrum of 
the spin waves: 

0(g) = y/(l + q‘^) (1 -h - cr^). (39) 

This dispersion relation is reproduced by the numerical 
simulations with a high accuracy, see Fig. 5(b). The crit¬ 
ical dependence, where the gap of the spectrum vanishes, 
reads 

Xc = \/l - (t2, (40) 

see thick dashed curve in Fig. 5(b). In the region between 
solid and dashed curves the normal state is metastable. 

The dispersion law (39) is symmetric with respect to 
the direction of the wave propagation: 0(g) = 0(—g). 
Unlike the easy-tangential case there is no effective mag¬ 
netic field A, because the curvature induced effective 
Dzyaloshinskii interaction is absent in the harmonic ap¬ 
proximation, cf. (28). The reason is that the equilibrium 
state is magnetized exactly in the normal direction Cn, 
which causes the degeneracy with respect to the sign of g. 
A similar behavior is known for thin films in the presence 
of Dzyaloshinskii-Moriya interaction, where the asymme¬ 
try in the spin wave spectrum vanishes if the system is 
saturated perpendicularly to the film plane. 


B. Easy—binormal anisotropy 

If the anisotropy axis is directed along Cb, one has the 
easy-binormal anisotropy, see (9b). The magneti¬ 

zation of the homogeneous (in the curvilinear reference 




(a) Quasi-binormal state: x = a = 1.5 




(b) Onion state: x = 1, cr = 2 


FIG. 7. (Color online) Magnetization distribution in the helix 
wire with C = -hi and easy-binormal anisotropy according to 
simulations data, see Sec. VIA. 


frame) state reads 


tan 20^ 


2C>z:cr 


coscj)^ = C = ±1, (41) 


Explicitly 0^ reads 


0*' = I [1 + sgn (Ccr)] - arctan 

q- = 1 + + q, 

Vi = Y^(l + — cr^)^ + 4x^(7^. 


(42) 


The magnetization of this state is close to the direction 
of the helix axis, hence we name it quasi-binormal state, 
see Fig. 7(a). It corresponds to the hollow-bar magneti¬ 
zation distribution in the helix microcoils.For different 
magnetization distributions see also Table I. 

The energy of the axial state reads 




1 - 

2w‘^ 


(43) 


Let us mention the formal analogy between the energy 
the “potentials” Vq , for the quasi-binormal state 
and the corresponding expressions [cf. (19)], Vq, U/ 
[cf. (17)] for the quasi-tangential state: the expressions 
for the quasi-tangential state can be used for the quasi¬ 
binormal one under the replacement k ^ a. 

The analogy between two states becomes deeper if we 
use another parametrization for the magnetization m 


m = cos 0 Ct — sin 0 sin 4> Cn + sin 0 cos 4> , (44) 


where 0 = 0(s) and 4> = ^>( 5 ) are the angles in the 
Frenet-Serret frame of reference: the polar angle 0 de¬ 
scribes the deviation of magnetization from the tangen¬ 
tial curve direction, while the azimuthal angle 4> corre¬ 
sponds to the deviation from the binormal. Similar to 
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(9), one can rewrite the energy terms as follows (cf. Ap¬ 
pendix A from the Ref. 11 for details): 


^ex = [0^ — /^sin$]^ + [sin0($' + “ KcosQcos^f , 




EB 


sim 0 cos^ ^ 


(45) 


Now one can easily see that the energy functional of the 
easy-tangential magnet transforms to the energy func¬ 
tional of the easy-binormal magnet under the following 
conjugations: ^ ^ 0, 0 ^ 4>, and >c ^ a. 

Similarly to the easy-tangential case, there exist 
two equilibrium states: the homogeneous state (quasi¬ 
binormal) and the periodic onion solution, see Fig. 7(b). 
The phase diagram, which separates these two states, is 
plotted in the Fig. 4(c). 

Now we discuss the magnons for the easy-binormal 
case. In analogy with the easy-tangential case, the 
linearized equations can be reduced to the generalized 
Scrodinger equation (27a) with the following “poten¬ 
tials” : 

,_l + x^+a^ + V,^ 

^2 — 2 ’ 

(46) 

A’’ = -KCOS0’’ - aGsind’’ = -KGVh 


where c^o = dTryM^, a is the damping coefficient, £ is the 
dimensionless energy, normalized by AttM^As^ with As 
being the sampling step of the natural parameter s. We 
consider four contributions to the energy of the system: 

The first term in Eq. (49a) is the exchange energy 

^ (49b) 

with £ = yjA/ (47rM|) being the exchange length. The 
second term determines the uniaxial anisotropy contri¬ 
bution 


n ^ 

£"" = -|E("^* •<")', (49c) 

where is the coordinate dependent unit vector along 
the anisotropy axis and Q is the quality factor, see (4). 
The third term determines interaction with the external 
magnetic field b 


N 

E^=-'^bi-mi, (49d) 

i=l 


The dispersion law has formally the form (30) with the 
corresponding “potentials” described above. The disper¬ 
sion curve is plotted in the Fig. 5(c) for some typical 
parameters, it is confirmed by the numerical simulations. 
The critical curve Xc(cr), where the gap of the spectrum 
vanishes, can be found numerically using condition (32). 
The critical curve >Cc{(j), calculated numerically is plot¬ 
ted in Fig. 4(c) (dashed red curve). For the approximate 
description of the critical dependence we use the trial 
function 



which fits the numerical results of Fig. 4(c) with an accu¬ 
racy of about 2 X 10“^, see the dotted curve in Fig. 4(c). 
In the region between solid and dashed curves the quasi¬ 
binormal state is metastable. 


VI. SIMULATIONS 

In order to verify our analytical results we numerically 
simulate the magnetization dynamics of a helix-shaped 
chain of discrete magnetic moments with i = 1, V. 
The form of the chain is described by Eq. (11b). Mag¬ 
netization dynamics of this system is determined by the 
set of Landau-Lifshitz equations 


where bi is the dimensionless external field, normalized 
by 47rMs. 

The last term in (49a) determines the dipolar interac¬ 
tion 

pd _ (As)^ -^' nii ■ jUj _ ^ {jUi ■ Vij) {rrij ■ rjj) 

Stt |r,,|3 |n,|5 • ^ ^ 

where Vij = 7^ - 7j. 

The dynamical problem is considered as a set of 3N 
ordinary differential equations (48) with respect to 3N 
unknown functions '^i(^) with i = 1, V. 

For a given initial conditions the set (48) is integrated 
numerically. During the integration process the condition 
\mi{t) \ = 1 is controlled. 

We considered the helix wire with length L = 500As, 
the exchange length £ = 3As and the quality factor Q = 2 
are fixed. The curvature k and the torsion r were varied 
under the restriction kAs/tt 1. 

In most of simulations we neglect magnetic dipolar in¬ 
teraction and consider the Heisenberg magnet with the 
energy 

£H ^ gex gan 

A. Equilibrium magnetization states 


1 drrii 

ujQ dt 


= rrii X 


98 

--h arrii x 

orrii 


rrii X 


98 ~ 

dmi_ 


(48) 


We start our simulations with easy-tangential mag¬ 
nets. In Sec. Ill A we found that the curvature and the 
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torsion causes the deviation of the magnetization from 
the anisotropy direction, which results in the magnetiza¬ 
tion distribution (17); such results are presented in Fig. 2 
by the curves for three different value of the reduced cur¬ 
vature K = 0.1,Xc, 1 in the wide range of the torsion 
G G (0;2). In order to verify our theoretical predictions 
we simulate numerically Landau-Lifshitz equations (48) 
in overdamped regime {a = 0.1) during a long time in¬ 
terval At ^ (ofCJo)”^. 

Numerically we model the anisotropic Heisenberg mag¬ 
net with the energy (50) and Q = 2. Simulation data are 
presented in Fig. 2 by filled symbols and labeled as (a); 
one can see an excellent agreement between out theory 
and simulations. The typical magnetization distribution 
is shown in Figs. 3(a), (b) and (c) for the quasi-tangential 
states and in Figs. 3(d) for the onion state. 

We also perform simulations for other anisotropy 
types. The magnetization distribution for the helix wire 
with easy-normal anisotropy is presented in Fig. 6(a) for 
the normal state and Fig. 6(b) for the onion one. For 
the case of easy-binormal anisotropy one has two possi¬ 
ble states: the quasi-binormal one [see Fig. 7(a)] and the 
onion one [see Fig. 7(b)]. 

The second stage of our simulations is to find the equi¬ 
librium magnetization state of a given helix wire. Numer¬ 
ically we simulate Eqs. (48) as described above for five 
different initial states, namely the tangential, onion, nor¬ 
mal, binormal, and the random states. The final static 
state with the lowest energy is considered to be the equi¬ 
librium magnetization state. We obtain that for each 
type of anisotropy the equilibrium state is either onion 
one or anisotropy-aligned state (quasi-tangential, nor¬ 
mal and quasi-binormal state for easy-tangential, easy- 
normal and easy-binormal anisotropy, respectively). We 
present simulations data in Fig. 4 by symbols together 
with theoretical results (plotted by lines). One can see a 
very good agreement between simulations and analytics. 


B. Dispersion relations 

For each anisotropy-aligned equilibrium state the 
magnon dispersion relation is obtained numerically. It 
is carried out in two steps. In the first step the helix wire 
is relaxed in external spatially nonuniform weak magnetic 
field 


bl = boef cos Sik^ 

for a range of wave-vectors = j/(300As) with j = 
0, 300. Here bo 1 is the field amplitude, Si = {i — l)As 
is position of the magnetic moment m^. The coordinate 
dependent unit vector ef determines the magnetic field 
direction: ef = for the quasi-tangential state and 
ef = Ct for normal and quasi-binormal states. 

In the second step we switch off the magnetic field and 
simulate the magnetization dynamics with the damping 
value a = 0.01 close to natural one. Then the space-time 


Fourier transform is performed for one of the magneti¬ 
zation components (we consider normal component for 
the quasi-tangential state and tangential component for 
other two equilibrium states). The frequency O which 
corresponds to the maximum of the Fourier signal is 
marked by a symbol for a given wave-vector = wk^ ^ 
see the top row of the Fig. 5. The absence of additional 
peaks in the spectrum is demonstrated by the dispersion 
maps below, see bottom raw of Fig. 5. 


VII. DISCUSSION 


We have performed a detailed study of statics and lin¬ 
ear dynamics of magnetization in the helix wire with dif¬ 
ferent anisotropy. We have limited our study by hard 
magnets, which can be well described by the model of 
anisotropic Heisenberg magnets. Our study was limited 
by the condition (5). 

Let us discuss how our model can be generalized 
with account of the long range magnetostatics effects. 
The non-local magnetostatic interaction for thin wires 
of circular and square cross-sections is known^^ to be 
completely reduced to a local effective easy-tangential 
anisotropy. It is important that such a conclusion sur¬ 
vives for the case of curved wires.Thus the magneto¬ 
static interaction can be taken into account as additional 
anisotropy. In general, one has to consider the model 
of biaxial magnet. Here we limit ourselves by the helix 
wire with easy-tangential magnetocrystalline anisotropy. 
In this case the magnetostatics effects can be taken into 
account by a simple redefinition of the anisotropy con¬ 
stants, leading to a new magnetic length 


K ^K^^ = K + ttM^, 


gAF 

w 



2i 

\/I -h 2Q 


(51) 


Thus our model (6) is also suitable for thin wires made 
of a magnetically soft material under the restriction 

h i, -, -. (52) 

K, T 

In order to check our predictions about effective 
anisotropy we perform numerical simulations with ac¬ 
count of the nonlocal dipolar interaction as described in 
Sec. VI. Numerically we integrate Eqs. (48) with the en¬ 
ergy (49). 

Eirst, we simulate the anisotropic wire with account 
the dipolar interaction with the energy (49). In this case 
we need to modify the magnetic length according to (51). 
Thus we also need to redefine the reduced curvature and 
torsion as follows 






(53) 


Eor the case (5 = 2, one get i = w and = 2^/\/5. One 
can see that we have a very nice agreement between the 
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analytical results (17) and simulations data, see yellow 
symbols in Fig. 2; we label these data as (am). 

The second kind of simulations with account of the 
dipolar interaction was aimed to verify the validity of 
our approach for soft magnets with Q = 0. For this pur¬ 
pose we model the soft isotropic wire with account of the 
dipolar interaction. According to (51) we get = 2£. 
Simulations data are presented in Fig. 2 by dotted sym¬ 
bols [labeled as (m)] for curvature and torsion redefined 
according to (53). By comparing simulations data with 
analytical results one can see the pretty good agreement 
in the wide range curvatures and torsions. Our simu¬ 
lations data for soft magnet differ from the theoretical 
predictions for hard magnets only for relatively high cur¬ 
vature in the vicinity of the boundary with the onion 
state. 

Thus we can conclude that our model of anisotropic 
Heisenberg magnet is physically sound also for thin wires 
made of a magnetically soft material. 

In conclusion, we have presented a detailed study of 
statics and linear dynamics of magnetization in the he¬ 
lix wire. We have described equilibrium magnetization 
states for three types of uniaxial anisotropy, according to 
possible curvilinear directions. All three cases have been 
realized experimentally in rolled-up ferromagnetic micro¬ 
helix coils.We have calculated the phase diagram of 
possible states in case of easy-tangential anisotropy: the 
quasi-tangential configuration (17) is energetically prefer¬ 
able for the strong anisotropy case. In this case the de¬ 
viations from the strictly tangential direction (corkscrew 
orientation^^) are caused by the torsion, the direction 
of the deviation depends on both helix chirality and 
the magneto chirality of the magnetization structure, see 
Eq. (18). In case of high curvature there is the onion equi¬ 
librium state (20) in analogues to the onion state in mag¬ 
netic ring wires^^’^^. The magnetization distribution (41) 
of the quasi-binormal state is directed almost along the 
binormal (hollow-bar orientation^^). In contrast to the 
quasi-tangential state and quasi-binormal one (which 
are realized for the easy-tangential and easy-binormal 
magnets, respectively), the normal state for the easy- 
normal magnets has several peculiarities: (i) it has the 
form of exact normal magnetization distribution along 
the normal direction Cn, see (35); (ii) the normal state 
phase is realized for small curvatures and torsions only: 

+ cr^/cTo < 1, see Fig. 4(b); (hi) the spectrum of 
spin waves on the normal state background is symmetric 
with respect to the direction of the wave propagation. 

The torsion of the wire manifests itself in the magne¬ 
tization dynamics: an effective magnetic field, induced 
by the torsion breaks the mirror symmetry with the spin 
wave direction. The dispersion law of spin waves (30) is 
essentially affected by this field. 

There is a connection between the helix geometry and 
the tube one: when the helix pitch vanishes, we have a 
close-coiled solenoid magnet, which properties are simi¬ 
lar to the thin shell nanotube. The spin-wave spectrum 
in the nanotube is known^^ to have a gap, caused by the 


curvature. This conclusion is in agreement with the dis¬ 
persion law for the helix wire, see Fig. 5(a). One has to 
note that the analogy between two systems is adequate 
under the restriction of vanishing torsions (cr ^ 0); this 
explains the absence of the linear shift in the dispersion 
law for the nanotube in comparison with (30). In gen¬ 
eral the transition from ID systems to 2D requires more 
accurate account of the dipolar interaction. 

We considered the simplest example of the curved wire 
with constant curvature and torsion. Our results can 
be generalized for the case of variables parameters hz{s) 
and t(s). To summarize we can formulate few general 
remarks about the curvature and torsion effects in the 
spin wave dynamics. The linear magnetization dynamics 
can be described by the generalized Scrodinger equation 
(27). In case of the straight wire, one has the standard 
Scrodinger equation for the complex magnon amplitude 

with the typical potential scattering. Th curvature 
induces an additional effective potential, the ‘geometri¬ 
cal potential’.This is described by the modification of 
effective potential U in Eq. (27b). Besides, there is a 
curvature induced coupling potential W : the problem be¬ 
comes different in principle from the usual set of coupled 
Scrodinger equations, see the discussion in Ref.^^. Due to 
the torsion influence there appears an effective magnetic 
field. The vector potential of this field is constant for 
the helix wire, see (26), hence the effective magnetic flux 
density B = 'V x A vanishes. Nevertheless the presence 
of magnetic field with the vector potential A breaks the 
mirror symmetry of the problem: the motion of magnetic 
excitations in different spatial direction is not identical. 

Let us mention the connection between the vector po¬ 
tential and the effective Dzyaloshinskii interaction: the 
total energy of the Dzyaloshinskii interaction oc 

f dsA-j with the current j = arg?/;, see Eq. (28). 

Using an explicit form of the integrand one can find that 
^ex ^ which reflects the relation between the topol¬ 
ogy of the wire (namely, helix chirality) with the topology 
of the magnetic structure (namely, the magneto chiral¬ 
ity). In this context it is instructive to note that there is 
a deep analogy between the Dzyaloshinskii-Moriya inter¬ 
action and the Berry phase theory^^. 

We expect that our approach can be easily general¬ 
ized for the arbitrary curved wires, where all potentials 
becomes spatially dependent: U(s), W{s), and A{s). De¬ 
pending on the curvature and the torsion these potentials 
can repel or attract magnons. In latter case there can ap¬ 
pear a well with possible bound states, i.e. local modes. 
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Appendix A: Onion-state solution 

We start from the static form of Landau-Lifshitz equa¬ 
tions (12): 

F{e,4>) = o, G{e,4>) = o (ai) 

with F and G being the nonlinear operators, 


G as follows 

N 

= '&n-, <y?l, • • • , ^n) COs(2n - l)x, 

n=l 

N 

G{0,(t)) = 'FGn (t?i,...,'(9„;<y3i,...,<y?„)sin2nx. 

n=l 

(A4) 

Here Fn and Gn are polynomials of the order N with re¬ 
spect to 'd/c and (/?/.. Then the Landau-Lifshitz equations 
(Al) results in the set of nonlinear polynomial equations 


F (6>, 0) = —d^^O — a cos 0 (x cos 20 — 2d^(p sin^ 0) 
-\-smOcosO (x + 5^0)^—(1-1-cr^)cos^ 0 , 

G {0, (j)) = sin^6> [~^xx^ + + cr^)sin (j) cos (j) 

—2ad^0 cos (j)] + sin 0 cos 0 [hcf sin (j) — 2d^0 {k + d^cj))] . 

(A2) 

By substituting here the expansion (20) in the form 


N 


^ix) = cos(2n - l)x, 


n=l 

N 


(A3) 


Hx) = -X + sin 2nx, 


n=l 


and expanding results into series over e up to the A-th 
order, one get the Fourier expansion of operators F and 


Gn {"Oi ^ ^ "0(pi ^ ^ (Pn) — 0 , 


n = 1, A, (A5) 


which can be solved numerically on lOk and pk with any 
precision. 

In order to calculate the energy of the onion state, we 
substitute the magnetization angles 0 and (j) in the form 
(A3) into the energy density (9), expand the results over 
5 up to the 2A-th order and average the result over the 
helix period. 




27r 


(A6) 


^ — ^ex + ^an ~ ^ ('^ 1 ? • • • 5 '^n; • • • 5 ^n) • 
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